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We use an extended two-band Kondo lattice model (KLM) to investigate the occurrence of differ- 
ent ( ant i-) ferromagnetic phases or phase separation depending on several model parameters. With 
regard to CMR-materials like the manganites we have added a Jahn- Teller term, direct antiferromag- 
netic coupling and Coulomb interaction to the KLM. The electronic properties are self-consistently 
calculated in an interpolating self-energy approach with no restriction to classical spins and going 
beyond mean-field treatments. Further on we do not have to limit the Hund's coupling to low or 
infinite values. Zero-temperature phase diagrams are presented for large parameter intervals. There 
are strong influences of the type of Coulomb interaction (intraband, interband) and of the important 
parameters (Hund's coupling, direct antiferromagnetic exchange, Jahn- Teller distortion), especially 
at intermediate couplings. 



I. INTRODUCTION 



The ferromagnetic Kondo lattice model, also known as 
double exchange or sd-model, is one of the basic mod- 
els in solid state physics. It is valid for systems which 
can be divided into two subsystems. The first one de- 
scribes the itinerant electrons and the other one local- 
ized electrons giving rise to finite permanent magnetic 
moments. A very prominent class of such materials are 
the manganites like Lai_a^Sra^Mn03 or Lai_a;Caa;Mn03. 
Due to a crystal field splitting the five 3d spin-up or- 
bitals of the manganese ion are split into three t2g and 
two Cg orbitals. The t2g spin-up states are fully occu- 
pied and provide the localized spin S = ^ and in the Cg 
states are n = 1 — x itinerant electrons. But the KLM 
alone does not explain the complex phase diagrams of 
the manganites^^'^ or other effects like the colossal mag- 
netoresistance (CMR)^. Thus other physical effects seem 
to be important. This can be the superexchange which 
takes place between the t2g electrons and leads to a direct 
antiferromagnetic coupling Jaf- According to the Jahn- 
Teller theorem a crystal with degenerated states breaks 
its symmetry and therefore lowers the energy. This re- 
sults in a splitting of the Cg orbitals. Last but not least 
the electrons in the narrow Cg bands experience a large 
Coulomb repulsion. 

The large number of competing interactions leads to ver- 
satile impacts on the phase diagran^. Not only that 
each effect influences the magnetic order for itself but 
they also act on each other. Thus it is necessary to un- 
derstand these interactions and their mutual effects in 
detail. Theoretical work on this has already been done 
using different met ho d^^^^ i ^ Q i ^ ^ . 

The composition of this paper is as follows: In the next 
section we will present our complete model. The meth- 
ods we have used to solve this model are described in 
Sec. mil Afterwards we will show the numerical results 
in Sec. [iVl There the magnetic phase diagrams will be 
presented which are derived by a comparison of the inter- 
nal energies of the different phases. Finally we will come 



to a conclusion in Sec. |Vl 

II. MODEL 

We use a two-band ferromagnetic Kondo lattice model 
as the main part of our Hamiltonian. 

HkLM = ^ TijC'l^^Cjaa " Jh ^ CTia ' (1) 

The first term describes the next-neighbor hopping of 
the itinerant electrons with the hopping matrix Tij and 

the fermion annihilation (creation) operators c^^^ for 
electrons with spin a in the band a . An on-site cou- 
pling of the spin of the conduction electrons (cr^a) to 
the spins of the localized electrons (S^) is done via the 
Hund's coupling Jh- The KLM for itself has a rich phase 
diagram^^'^^ and is often used as the only part to de- 
scribe real material^i^^i^ with some success. But spe- 
cial additions are needed to explain special features or 
more complex materials^. One of the most prominent 
interactions is the Coulomb repulsion represented by the 
Hubbard part 

Hu = riiaariiaa + ^ U"^"^ niaaTliaa' , (2) 

where riiaa = ctaa^ioio- and a bar above an index means 
the opposite band or spin. If we choose U^^ = there is 
only mtraband repulsion and with U^^ 7^ there is also 
interband repulsion. Besides this direct electron-electron 
interaction there can be a direct coupling of the localized 
spins. 

Haf = Jaf ^ Si • Sj (3) 

This interaction is often generated by the superexchange 
of electrons. Therefore it is only used as an antiferromag- 
netic coupling {Jaf > 0) in this paper. The last exten- 
sion is due to the Jahn- Teller effect (JTE), which lowers 
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the degeneracy of electron states by reducing the symme- 
try of the crystal. Mostly it is related to Sd-orbitals, e.g. 
the e^-orbitals in the manganites. In a standard notation 
it is written as 



Ql)- (4) 



The Q2(3)i are special JT modes and the T- are pseu- 
dospin operators where the band index replaces the spin 
index of the usual spin operators. In our paper we set 
/cjT = 1. With this final part we have the complete 
Hamiltonian 



W = ^KLM + Hu + Haf + HjT 



III. CALCULATION METHODS 



(5) 



To get magnetic phase diagrams at zero temperature 
we have to calculate the internal energy of the differently 
ordered magnetic configurations. We primarily focus on 
the ferromagnetic alignment and get the antiferromag- 
netic phases by a division of the whole chemical lattice 
into ferromagnetic sublattices. 



A. Ferromagnetic Phase 

For the internal energy we need the quasi-particle den- 
sity of states paa{E) (QDOS) which we can get from the 
one-particle Green's functions 
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In principle these Green's functions can be specified by 
solving the according equations of motion (EOM). Unfor- 
tunately there is no known exact analytical solution for 
this model thus we need approximation methods. Just 
looking at the first two parts in ([5]) we can directly use 
the interpolating self-energy approach (ISA)^^ which has 
been successfully applied for the description of real mate- 
rials (e.g. Ref.^^^^^^^^). This approach fulfills the exactly 
solvable limiting cases of the KLM (ferromagnetically or- 
dered semiconductor, atomic limit, second order pertur- 
bation theory) and interpolates them by fitting free pa- 
rameters via a high energy expansion. Therefore we ex- 
pect reasonable results even between the limiting cases 
and it should hold for all band occupations, temperatures 
and all orders of Hund's coupling. Within this approach 
we get the self-energy 
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containing 
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In addition to the model parameters of ([5]) we need A^a 
and {riaa) which can be self-consistently calculated via 
the spectral theorem from the full Green's function 



^ + /i-T,.(k)-SjSA(^) 
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with the Fermi function f-{E). The expectation value 
(S^) has to be considered as an external parameter or can 
be taken from another method^^. In the upper Green's 
function the Hubbard part of the Hamiltonian is incorpo- 
rated in the j^a via an effectice medium approach. Nor- 
mally the Green's function ([8]) contains a second part 
(upper Hubbard band) which we left out choosing re- 
spectively > J or C7^^' > W, J. Thus (O only de- 
scribes electrons which have no Coulomb repulsion part- 
ner at their site. The probability that an electron has no 
repulsion partner is 



1 - {ria-a) -{n-a,a) - {ri-a-a) • 



(11) 



Since the jaa are in the numerator of (|8]) they infiuence 
the spectral weight of the Green's function but they also 
act on the bandwidth in 

T„.(k) = T(°) + 7.. (e(k) - Ti°)) (12) 

where T^^ are the centers of gravity of the bands. For 
this reason we have a total spectral weight 



4 — n intraband 

4 — 3n intra+interband 



(13) 



which depends on the electron density n. Especially for 
n = 1 and interband repulsion we have completely filled 
bands (F = n = 1) and therefore a Mott insulator. 
We now have to add the Jahn- Teller Hamiltonian of 
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Eq. This can be easily done when we treat 

the phonon operators classicahy and use a mean-field 
approximation^' After performing the mean-field de- 
coupling the phononic variables Q2i-, Qsi are only coupled 
to mean values of the electrons and the ground state is de- 
fined by the mean values (Q2i), {Qsi) due to the absence 
of quantum fiuctuations. Introducing spherical coordi- 
nates Q2 = Qcos^, Qs = sin^ and dressed operators 



^iaa = — l — ^' ^ (cos Cj^a,3z'^ — r 



.2 -fsm-Q^^^2. 



Ciaa=+1 = e*^^^(sin-Q^,322_^2 - COS-Q^,3,2_^2Xl4) 

we can replace the two modes in (|4]). These operators 
are superpositions of 3(i-orbitals in z-direction or in the 
x-?/-plane, respectively. If we also assume translational 
invariance (i.e. a non-cooperative JTE) the JTE is not 
dependent on the angle but only on the magnitude of 
the distortion Q. In the ground state this distortion is 
defined by Q = g{An) and we get the new form of the 
Hamiltonian (|4]) 



ny , (15) 
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Figure 1: (color online) The different types of magnetic order 
which are investigated in this paper. Arrows show the spin 
direction. A,C,G-type antiferromagnetism, FM/PM ferro- 
/ paramagnetism. 



the Hamiltonian ([5]) are local thus we have to sum over 
one more index u which distinguishes the two sublattices. 
The hopping part needs two indices v' since there is of 
course a hopping between different sublattices. Neglect- 
ing again Haf we get 



lA 5 



(17) 



where H^^ includes all electronic interaction parts (dou- 
ble exchange, Hubbard and JTE) in the magnetic sub- 
lattice v and 



vv' + 
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with 2;ce=±i = ±1 and the occupation difference 
(An) = ^ ((nc=-i,a) - 



(16) 



which has to be calculated self-consistently using (|9]). 
Since we want to add this term to the Green's function 
([8]) we can neglect the second term in ([15]) because it 
plays no role in the equation of motion. But we have to 
keep it in mind if we calculate the internal energy. The 
first term means only a shift of the energy of the two band 
in different directions. Hence it can be incorporated in 
the centers of gravity in (p!2|) which have then to be 
calculated self-consistently, too. 

Now we have found a solution for all parts of the Hamil- 
tonian ([5|) which contains electron operators and thus we 
can calculate the quasi-particle DOS p^cr. The treatment 
of the superexchange part (|3]) will be described in IIII CI 



The dispersions e 
integrals 



are the Fourier-transformed hopping 



ik(Rr-R,^ ) 



(19) 



We only allow hopping to next neighbors in the chemical 
lattice. For solving the equation of motion (EOM) we 
define the Green's function 



Gi::;:(i?) = ((ck„..;c+„,,,)) . 



(20) 



To get the quasi-particle DOS p^^^iE) we have to calcu- 



late the EOM for the G{^^^(-B) which is 
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B. Antiferromagnetic Phases 

In the section above we found a solution for the fer- 
romagnetically ordered system. To extend this result 
to antiferromagnetic order we divide the whole chemical 
lattice into ferromagnetic sublattices. This means that 
we assume a Neel state which is known not to be the 
ground state due to quantum fiuctuations. But at least 
for spins 5 > | it should be close enough to the ground 
state to make reliable conclusions^^^^^. For simplicity we 
only investigate types of antiferromagnetism which can 
be divided from a simple cubic lattice into two sublat- 
tices (Fig. [1]). All of the electronic interaction parts in 



Here [. . . ] represents the commutator, v means the 
opposite sublattice. The higher Green's function is 
only affected by local interactions therefore we can 
use the sublattice self-energy ^c\^a.cjv,ii\p\\c^^^J) = 
Ml^aAE)GZa{E)^'' Thus M{^^^{E) is the self-energy 
of a pure ferromagnetic case comparable to (|8]). Now 
we need the Green's function G^^^{E) coming from the 
interlattice hopping. This leads to a second EOM 
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Ml^A (23) 



Other simplifications can be done if one considers some 
symmetries. The intra-sublattice hoppings should be the 
same for both sublattices because of the same chemical 
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structure, i.e. = = ek- This also holds for the 
inter-sublattice hoppings = = t\^. Since the two 
sublattices only differ by the spin direction we can replace 
the self-energy of the opposite sublattice by switching the 
spin 



(24) 



For the formal self-energy Mk^a we now use the ISA self- 
energy plus the terms describing the Jahn- Teller effect 
and the Coulomb repulsion like in Eq. ([8]) . With ([22|) and 
([23|) we find the full Green's function of one ferromagnetic 
sublattice 
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{E)) 



(25) 



where T«,(k) = T^^^ + 7«,(gk - TD similar to 
The Hubbard interaction is contained in the and the 
JT splitting in the T^^^ = Za{An)g'^ (cf. Eq. ([15])). 
Although this treatment needed some approximations it 
fulfills the limiting cases Jh and (S^) ^ 0. In these 
cases the paramagnetic DOS is reproduced. 



(0)n 



C. Internal Energy 

To decide which magnetic phase is preferred at a spe- 
cial parameter set we have to compare the internal ener- 
gies. In the KLM it is simply given as 



_ {Hklm) 
N ~ N 



+ CO 



(26) 



It is a special feature of the KLM that some parts can- 
cel each other out and the typical e(k)-part is missing. 
The electronic parts of the Hamiltonian ([5]) which do not 
belong to i^KLM are contained in the Green's functions 
dH]) and ([25]) via an effective medium approach. Thus 
they are already taken into account by the quasi-particle 
DOS paa{E). We just have to add the second term of 
(p!5|) which had no influence on the EOMs but on the 
internal energy. This means for the electronic internal 
energy 



Uei/N = U'/N+^g^Anf 



(27) 



The only missing part is the energy of direct coupling 
between the localized moments (j3j). Since we are in the 
zero-temperature regime we use a mean-field decoupling 

-^""^^^ =JAF{S^f{a,-a^,). (28) 



N 



N 



We defined a(a)p as the number of (anti-)parallel aligned 
next-neighbor spins. This energy is lowest for highest an- 
tiferromagnetic order (G-type, all neighbors are antipar- 
allel aligned) and vice versa highest for ferromagnetic 



order. The whole internal energy is now 
U/N = Uei/N + Uaf/N 



(29) 



Ferromagnetic saturation is assumed for the respec- 
tive (anti-) ferromagnetic phases. Thus we set (S^) to 
the maximum value which is in the ISA^^ (S^)^^^ = 

\ ' IIlclX 

S^^^Y^' For the paramagnetic phase there is a vanishing 
mean value of the magnetization {S^) = and therefore 
also no contribution to the spin's internal energy (E 



D. Phase Separation 

In the regions of phase separation the complete internal 
energy can be expressed by the internal energies of the 
phases A, B as 



U 
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{l-y)U^{nA)^yU''{nB) 



(30) 



with the volume fraction y covered by phase B. Each 
phase has its own electron density nA,B and (1 — y)nA + 



Ut 




n 



Figure 2: Schematic view on the Maxwell construction. Be- 
tween ni and n2 the internal energy curves of the phases A 
and B are replaced by a straight line representing the occur- 
rence of phase separation. 
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yuB = n. Minimization of U^^ according to ua^b leads 
to conditions for the boundaries of the phase separated 



area 



driA 



dni 



(31) 



Between ni and n2 the phases A and B coexist (cf. Fig. 

ED- 



IV. RESULTS 

In this section we wih present the numerical results of 
the self-consistent calculations. The effect of the different 
extensions will be studied in detail. All results are given 
for the spin 5* = |. 



A. Coulomb Interaction 

Basically the Coulomb repulsion can act in three ways. 
It can be turned off completely, act only between elec- 
trons in the same band or between electrons of differ- 
ent bands, too. In our treatment it affects the band- 
width and the spectral weight of each band. The orig- 
inal values are reduced by the factor ^o,a- in (HH)- Be- 
cause these factors are dependent on the occupation num- 
ber the main differences occur at high electron densi- 
ties. For example at n = 1 the bands are quarter-filled 
if there is no Coulomb repulsion, one-third-filled at in- 
traband repulsion and completely filled with additional 
interband repulsion. Ferromagnetism at large couplings 
Jh in the pure double-exchange model is most favored 
at quarter-fillings^^. This is the maximum filling that 
can be reached with an absent Coulomb interaction for 
densities < n < 1. On the other hand there is no fer- 
romagnetism in the DE model at half-filling where the 
chemical potential lies between the bands. A compara- 
ble situation in our extended KLM can only be achieved 
with intra+interband Coulomb repulsion dX n ^ 0.94 for 
(anti) ferromagnetic phases (Fig. [3|). So one can explain 
basic differences of the phase diagrams in Fig. [3l 
Without Coulomb interaction we have effective quarter- 
filling at n = 1. This means that ferromagnetism is 
mostly provided by the itinerant electrons and the en- 
ergy difference /SUei = U^^^ — U^i^ has its maximum 
at this electron density. With a finite antiferromagnetic 
coupling Jaf antiferromagnetic phases begin to appear 
at lower electron densities n where the AUei are small 
(cf. Fig. I8l9p . Thus the total energy differences are 
mainly governed by the energy of the localized spins ([28|) 
which is density-independent. The G-type AFM gains 
most energy from the direct antiferromagnetic coupling 
and exists at the lowest electron densities (compare Sec- 
tion [IVC]). 

If we have intraband repulsion antiferromagnetic phases 
also appear at higher densities. Because of the reduction 



of the total spectral weight we cross the effective quarter- 
filling and the absolute energy difference of the phases is 
reduced again (cf. Fig. [8|). Therefore we get antiferro- 
magnetic phases if we have a sufficient antiferromagnetic 
coupling Jaf at n = 1 without losing the ferromagnetic 
phase for n < 1. This was not possible without Coulomb 
interaction where the ferromagnetic phase stays longest 
at n = 1. 

When we also add the interband repulsion the spectral 
weight will be reduced even more. Thus the maximum 
absolute electronic energy difference for the ferromag- 
netic phase lies at n 0.4 (cf. Fig. [9|). There the 
ferromagnetic phase exists longest with increasing Jaf 
(cf. Fig. [3|). As a second feature the upper subband 
will be occupied for electron densities n ~ 0.94 (resp. 
n ^ 0.91 for the paramagnetic phase). Because of the 
extremely reduced spectral weight the energy differences 
become very small and we get many small regions of dif- 
ferent phases for n > 0.94 at Jaf = 0. This occurrence 
of many small (anti-) ferromagnetic regions could be a 
hint for phase separation or spin-canted states. We will 
discuss phase separation in section [TVE[ With finite Jaf 
the interval 0.94 < n < 1 is dominated by the energy of 
the localized spins and that is why the G-type AFM ap- 
pears. 

In general the Coulomb interaction favors ferromagnetic 
order. For both types of repulsion the whole phase di- 
agram would be ferromagnetic if Jaf = (except at 
n ^ 0.94 at interband repulsion) even for very low Jh 
and U (Ref.^^). It can be argued that the interplay be- 
tween the two ordering mechanisms, the Coulomb inter- 
action and the double exchange, makes it easier for the 
system to go to the symmetric state. But the reduction 
of the spectral weight and the bandwidth lower the ab- 
solute energy differences of the itinerant electron system. 
Thus one needs a smaller Jaf to create antiferromagnetic 
phases. For this reason the interband Coulomb phase di- 
agram in Fig. [3] has larger regions of antiferromagnetic 
phases compared to the others. 

From now on we will only discuss the cases of intra or 
intra+interband Coulomb repulsion. 



B. Hund's Coupling 

In the ferromagnetic KLM the Hund's coupling Jh 
in principle favors ferromagnetism. But because of the 
complex interplay of a parallel alignment of the local- 
ized/itinerant spins and the hopping there also occur an- 
tiferromagnetic phases. These exist mainly at low and in- 
termediate couplings Jh and at high electron densities^^. 
This difference between low and strong coupling has its 
origin in the splitting of the spin-up and spin-down bands 
(Fig. [3j). The quasi-particle DOS changes strongest at 
low couplings and there are more spin-down states occu- 
pied. At some Jf^ the main part of the spin-down band 
is shifted above the chemical potential depending on the 
electron density n. For strong couplings there is no major 
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Figure 3: (color online) upper line: QDOS and chemical potential (vertical lines) at n = 0.94 for ferromagnetic (solid black line) 
and A-type antiferromagnetic (dashed red line) order at Jh = 0.75eV. left: No Coulomb repulsion, middle: only intraband 
repulsion, right: intra+interband repulsion all: W = leV, Jaf = l.TmeV, g = 0. lower line: Phase diagrams in dependence of 
the Hund's coupUng Jh in eV and the electron density n for different types of magnetic order. For vanishing antiferromagnetic 
coupling Jaf the whole phase diagram for intraband repulsion would be ferromagnetic and also everything below n ^ 0.94 for 
intra+interband repulsion. 



variation of the shape of the QDOS. Only the spHtting 
gets larger with larger Jh- For example this is reflected 
in the Curie temperature which cannot be increased by 
increasing Jh at large couplings^^. 

Since we have an additional Hubbard part to our basic 
KLM we have two supporting mechanisms that create 
ferromagnetism. This means we would have an almost 
complete ferromagnetic phase diagram if we would choose 
Jaf = 0- With a finite Jaf one can see the emergence 
of different antiferromagnetic phases. Figure 2] shows the 
disparity between low and strong coupling. For low Jh 
the boundaries between the phases are dependent on the 
Hund's coupling and the electron density. On the other 
hand, at strong couplings, these boundaries occur at con- 
stant densities Uc and therefore they are almost vertical 
in the phase diagram. The phases are now very stable at 
certain densities n concerning a change of Jh- This sta- 
bility vanishes with increasing Jaf when the electronic 
energy differences become to small. 

The main effect of Jaf happens of course at small differ- 
ences of the electron energies (Fig. [9j). That is given, for 
intra+interband Coulomb interaction, at small and high 
n. A small n means that there are not enough electrons 
to create large absolute energy differences. On the other 
hand a large n leads to an occupation of the upper sub- 
band (cf. IIV Ap which is at higher energies. For n ^ 1 
the electronic internal energies for all phases go to zero 
and therefore the energy differences vanish, too. 



C. Direct Antiferromagnetic Coupling 

In the last sections we investigated the influence of elec- 
tronic correlations on the phase diagram. Their interplay 
mainly provided ferromagnetism. Thus we have already 
seen that we need a finite Jaf to get larger regions of an- 
tiferromagnetic phases. Compared to the former mecha- 
nisms the direct antiferromagnetic coupling has no influ- 
ence on the quasi-particle DOS. At finite temperatures it 
would of course act on the magnetization (Sz)- But as 
we look at zero-temperature behavior we have restricted 
(Sz) to its maximum value. That is why Jaf only influ- 
ences the energy of the localized spins (|28|) and can be 
examined more separately. 

The energy ([28|) is primarily dependent on the number 
of (anti)parallel ordered spins. As the G-type AFM has 
the most antiparallel neighbors it gains most from a fi- 
nite JaF' Vice versa the ferromagnetic order is most 
suppressed by it. The other phases get an intermedi- 
ate energy change depending on their structure (compare 
Fig. [1]). This energy change competes with the energy 
difference of the electronic subsystem. Thus the G-type 
emerges at low absolute energy differences AUei (at n 
or n ^ 1 for intra+interband Coulomb repulsion) even for 
low Jaf- The A- and C-type typically begin at interme- 
diate parameters Jaf and n where the absolute AUei is 
the largest (Fig. I5l9p . The G-type normally is unpre- 
ferred by the itinerant electron system in this region. 
It is not always possible to get all phases by varying Jaf • 
At some parameters sets only ferromagnetism and the G- 
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Figure 4: (color online) upper line: QDOS and chemical potential for different Hund's couplings Jh and n = 0.9. Solid (black) 
line ferromagnetism, dashed (red) line A-type AFM. lower line: Phase diagrams in dependence of the Hund's coupling Jh in 
eV and the electron density n for different types of magnetic order, left: Jaf = 0.56meV, middle: Jaf = 0.94meV, right: 
Jaf — l.SlmeV all: intra+interband Coulomb repulsion W — leV, ^ = 0. 



type appear, for example. Thus it seems that complex 
phases diagrams, like those of the manganites, need the 
interplay of different interactions. More phases can ap- 
pear if one has a finite electron-phonon coupling g (Fig. 
[5j) . As described in the next section this coupling has an 
unequal effect on the particular phases. The G-type can 
be reached everytime with a sufficient Jaf, of course. 



D. Jahn- Teller coupling 

The Jahn- Teller theorem says that a symmetric crys- 
tal with degenerated states lowers its symmetry so that it 
reduces its energy. Thus all the magnetic phases should 
lower their energy when the two orbit als (p!l|) split. But 
as long as we calculate the splitting (p!6|) self-consistently 
it is not clear which coupling Qc is necessary to lift the 
degeneracy. Figure [6] shows indeed that the critical cou- 
pling is not the same for the single ordering types. For 
example at low couplings ^ ^ 0.3 the G-type AFM ap- 
pears first in the density interval 0.6 < n < 0.9 when we 
use intra+interband Coulomb interaction. On the other 
hand the A-type AFM shows up first at large g and low 
densities. It is not only important which phase is JT- 
split but also how much a single phase profits from the 
splitting. The effect of the JT splitting on the electronic 
internal energies can be seen in Fig. [9l 
It was mentioned before that the phase diagrams with 
finite Coulomb interaction and vanishing Jaf are almost 
purely ferromagnetic. But with an increasing coupling 
g the G-type starts at n 0.8 and the crystal gets G- 



type ordered. When one increases g further on, the JT 
splitting starts in the paramagnetic phase and it is now 
preferred, even though the G-type is still split, too. With 
further increase of g other phases also appear. The en- 
ergy lowering due to the JT-part of the Hamiltonian (|4]) 
depends on the occupation difference (p!6|) . The maxi- 
mum value of (p!6|) is (An)^^^^ = n. Thus the JT energy 
is only sufficient to create AFM phase at higher densities 
(Jaf = 0). 

If we use a finite Jaf one sees that the single phases 
develop at the boundaries of their splitting zones with 
increasing Jaf- That is why the A-type AFM starts at 
low densities and large g where it is the only split phase. 
In the strong coupling region {2g'^n > W) the bound- 
aries are not electron density dependent any more. Thus 
these boundaries get vertical similar to the large Jh case 
in Sec. IIVBI But the edges are shifted compared to the 
non-split region at low g. 

This different behavior for either no, intermediate or 
strong couplings can also be seen in Fig. [5l At low and 
strong g we see no sharp edges in the phase diagrams 
which corresponds to slight changes in the electronic en- 
ergy ([27|) . This is because there is either no splitting or 
at large couplings the bands in all phases are split to their 
maximum value 2g'^n. For intermediate couplings there 
are again different critical densities Uc for each phase. 
That is why we get sharp edges at these values. 
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Figure 5: (color online) Phase diagrams in dependence of the direct antiferromagnetic coupling Jaf in eV and the electron 
density n for different types of magnetic order. Additionally there is either no, intermediate or strong electron-phonon coupling g. 
upper line: intraband Coulomb repulsion left: g — Ov^ey, middle: g — 0.9A/ey, right: g — l.O^/eV lower line: intra+interband 
Coulomb repulsion left: g — Ov^ey, middle: g = O.GA/eF, right: g — l.Ov^eF all: W = 3eV, Jh = 2eV. 



E. Phase Separation 

In contrast to the homogeneous phases presented in 
the sections above, it is also possible that we have phase 
separated regions in the phase diagram. To see where 
a mixture of two phases can exist we use the method 
described in section IIIIDi Indeed we find that there is 
always a phase separated area between two phases. As an 
example we show in Fig. [Tithe regions of phase separation 
for the three different types of Coulomb interaction. If 
we have no Coulomb interaction we can get a very broad 
range of phase separation, especially at large Jh and low 
n. With the increase of Coulomb repulsion partners the 
phase separation occurs only in smaller intervals around 
the original boundaries of the two phases. This has its 
reason in the gradient of the internal energy curves. If 
the gradient is small we get large regions of phases sep- 
aration due to the construction we made in Eq. (|3T]) . 
Figure [8] shows that we get the smallest gradient at low 
densities and the least influence of the Coulomb interac- 
tion. 

It is also possible that complete homogeneous phases can 
vanish due to phases separation. For example, the ferro- 
magnetic phase in Fig. [71 at intra+interband Coulomb re- 
pulsion is, except for a small stripe, covered by phase sep- 
aration of ferromagnetism/G-type antiferromagnetism or 
ferromagnetism/ A- type respectively. 
The inclusion of phase separation into the other phase 
diagrams leads to a qualitatively similar picture. Phase 



separation occurs always between two phases again in a 
more or less broad range. This would increase the va- 
riety of those phase diagrams even more and shows the 
complexity of the underlying model. 



V. SUMMARY 

We have investigated the influence of several exten- 
sions of the two-band Kondo lattice model on the zero- 
temperature magnetic phase diagrams. Those extensions 
were a Hubbard term, direct antiferromagnetic coupling 
of the localized spins and electron-phonon coupling due 
to a Jahn- Teller part. We used self- consistent calcula- 
tions, which go beyond the mean-field level. It has been 
seen that these extensions have a large impact on the 
phase diagrams. Because of the mutual effect on each 
other we can get very complex phase diagrams. 
The Coulomb interaction normally leads to a ferromag- 
netic state. But due to a lowering of the spectral weight 
and bandwidth of the quasi-particle density of states it 
also lowers the absolute energy differences of the elec- 
tronic subsystem. That is why the occurrence of anti- 
ferromagnetic phases needs a smaller antiferromagnetic 
coupling Jaf at intra+interband Coulomb interaction 
compared to the case that one has intraband repulsion, 
only. Especially because of the strong Coulomb interac- 
tion we used in our work we always needed a finite Jaf 
to create larger regions of antiferromagnetic phases. Big 
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Figure 6: (color online) Phase diagrams in dependence of the electron-phonon coupling ^ in veV and the electron density 
n for different types of magnetic order, upper left: Occurrence of the Jahn- Teller splitting. The lines mark the beginning 
((An) ^ |n) of the splitting, remaining three graphs: Favored magnetic phases for upper right: Jaf = (mind the different 
n-scale!),/oi(;er left: Jaf — 3meV and lower right: Jaf — 4.7meV all: Jh = 2eV, W = 3eV. The JTE lowers the energy. Thus 
new phases develop at the beginning of the JT splitting with increasing Jaf and lead to a complex picture for intermediate 
coupUngs g. Values for n > 0.94 are left out in the upper line due to graphical reasons. 




Figure 7: (color online) Phase diagrams of Fig. [3] but with the inclusion of phase separation (striped). The phase separation 
consists of the two phases which are adjacent to the left and right side of the phase separated region, left: no Coulomb 
interaction middle: intraband Coulomb repulsion right: intra+interband Coulomb. Zones of phase separation become smaller 
with the increase of the influence of the Coulomb interaction. 



values of Jaf always provide G-type AFM but interme- 
diate values can lead to other phases, too. 
The electron-phonon coupling g can have a very subtle 
influence on the magnetic ordering. We have to differ 
again between low, strong and especially intermediate 
couplings. At low g there is no splitting of the two bands 



and for large g the important band occupation difference 
([16]) is saturated at (An) = n. The intermediate cou- 
plings show a large variety of different phases. 
Further on we looked on the appearance of phase sep- 
aration. We found that there is indeed always a phase 
separated region between two homogeneous phases. De- 
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Figure 8: (color online) Differences of the electronic internal 
energy (|27p between ferromagnetic and paramagnetic phase 
for the single types of Coulomb interaction. The maximum 
of the absolute difference occurs at unequal electron densities 
n. At intra+interband repulsion the upper subband in the 
paramagnetic phase is occupied for n > 0.91 and for the fer- 
romagnetic phase at n > 0.94. Thus there are larger energy 
differences in this region except for n ^ \ where the bands 
are completely filled and the energy difference goes to zero. 
Parameters: Jh — 1.5eV, W — leV, g — ^ 



pending on the parameters this region can be very broad 
or very small. For special cases a originally homogeneous 
phase can be completely covered by phase separation. 
To get complex phase diagrams, like for the manganites, 
extension to the ferromagnetic Kondo lattice model seem 
to be necessary, especially at large Jh and Hubbard U . 
Here intermediate values of the direct antiferromagnetic 
and/or electron-phonon coupling play an important role 
to get a big diversity of magnetic phases. This could 
mean that approaches which are only valid in low or 
strong coupling regimes are maybe not appropriate to 
describe such material classes. 

The magnetic order can surely be influenced by other ef- 
fects. That could be orbital and charge order as well as 
the temperature, of course. These extensions would in- 
crease the complexity of the calculations much more and 
are left for later work. 



Appendix A: ELECTRONIC ENERGIES 



In this section the differences of electronic internal en- 
ergies (|27|) between the (anti) ferromagnetic phases and 
the paramagnetic phases are shown. It can be seen in 
Fig. [8] and [9] that the lowest absolute energy differences 
occur at low electron densities n. In our treatment of 
the internal energy of the local moment system ([28|) the 
whole internal energy would be equally shifted for all n 
at a finite Jaf- Thus it is most dominating at lower 
n for all types of Coulomb interaction. At intra+inter 
band Coulomb repulsion at n > 0.91 the paramagnetic 
phase is always unpreferred but the energy differences 
between the (anti) ferromagnetic phases are very small. 
Therefore the local moment internal energy is most im- 
portant here, too. The paramagnetic internal energy is 
not affected by a finite Jaf, the ferromagnetic phase and 
the A-type AFM are shifted to higher energies and the 
C and G-type AFM to lower ones. Therefore the G-type 
is preferred in those regions. In contrast to that we have 
a maximum of the absolute energy difference at a special 
^max- The position of this maximum is mainly depen- 
dent on the type of the Coulomb interaction. It is at 
n = 1 for the case of vanishing Hubbard repulsion and 
gets lower for finite Coulomb interactions. Near the nmax 
the ferromagnetic phase remains longest with increasing 
Jaf- 

If we have a finite JT coupling g the JT bands are able to 
split. Since this splitting is calculated self-consistently it 
occurs only for special electron densities depending on g 
and the type of the magnetic ordering (cf. Fig. [6]). As the 
JT splitting reduces the energy in most cases, magnetic 
phases with a finite band occupation difference (An) are 
preferred. Figure [9] shows the occurrence of the splitting 
for the single phases. Ferromagnetic order, for example, 
is only possible below n ^ 0.8 where the ferromagnetic 
system gets a finite (An) for the given parameters. At 
lower n the JTE is less important because the maximum 
value of (An)"^^^ = n is lower and the system is less 
profiting from a splitting. 
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